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Abstract. 

The electric resistance between two arbitrary nodes on any infinite lattice structure 
of resistors that is a periodic tiling of space is obtained. Our general approach is based 
on the lattice Green's function of the Laplacian matrix associated with the network. 
We present several non-trivial examples to show how efficient our method is. Deriving 
explicit resistance formulas it is shown that the Kagome, the diced and the decorated 
lattice can be mapped to the triangular and square lattice of resistors. Our work can 
be extended to the random walk problem or to electron dynamics in condensed matter 
physics. 
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1. Introduction 

In electric circuit theory one of the classic problems is to calculate the resistance between 
two arbitrary grid points on an infinite square lattice of identical resistors. This problem 
has been published by van der Pol and Bremmer [Tj. In a special case, it is well known 
that the resistance between two adjacent grid points of an infinite square lattice is half 
the resistance of a bond, and an elegant and elementary solution of the problem is given 
by Aitchison |2]. Over many years numerous authors have studied this problem and its 
extension to triangular, honeycomb and hypercube lattices (for some relevant references 
see, eg, ^W)- 

The resistor network can systematically be treated by the Laplacian operator of the 
difference equations governed by Ohm's and Kirchhoff's laws. Then the Green's function 
corresponding to the discrete Laplacian operator can be related to the resistance between 
two arbitrary nodes on a resistor network. This idea has been applied by Kirkpatrick 
to study percolating networks of resistors j6]. Recently, it has been demonstrated that 
the application of the Green's function is a very efficient way to find the resistance of 
an infinite resistor network as well [5]. The concept of a Green's function is widely used 
in the literature. Economou's book [7] gives an excellent introduction to the Green's 
function. Katsura et. al. have a review of the lattice Green's function [8] and for more 
works on this topic, see references in [5]. From the equation for the Green's function one 
can, in principle, derive some recurrence formulas for the resistances between arbitrary 
grid points of an infinite lattice (for square lattice of resistor network see reference [5]). 
The Green's function can also be a useful tool to study, eg, capacitor networks [9] or the 
resistance in a perturbed lattice in which one of the bonds is missing in the lattice [10] . 
The Green's function for the anisotropic diamond lattice was discussed by Hijjawi [H] 
and the analytical properties of the Green's function in all dimensions were investigated 
by Guttmann |12j . 

The problem of a finite resistor network is equally interesting in circuit theory. 
Recently, Wu has developed a theory to calculate the resistance between arbitrary nodes 
for a finite lattice of resistors ^13j, Tzeng and Wu for impedance networks [H]. The 
corner-to-corner resistance and its asymptotic expansion for free boundary conditions 
were obtained by Essam and Wu [I5] and for other boundary conditions the latter was 
calculated by Izmailian and Huang [T6] . 

Calculation of the resistance in electric circuit theory can also be relevant to many 
other problems, such as random walks [TTl [18] and first-passage processes [19j. This 
connection is based on the fact that the difference equations for the electrical potentials 
on a lattice point of a resistor network is the same as those that occur in the above 
mentioned problems. Thus the resistance problem can be regarded as a problem to 
solve the difference equations for an infinite network. 

In this work we generalize the Green's function approach developed in reference [5] 
for a resistor network that is a uniform tiling of d dimensional space with electrical 
resistors. In two dimensions, tiling is a collection of plane figures that fills the plane 
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with no overlaps and no gaps. Generalization to other dimensions is also possible. 
Figure [1] shows an example for a uniform tiling with resistors. Here we derive an 



Figure 1. An example for a lattice structure of infinite resistor networks discussed in 
this work. The lattice is a periodic tiling of the plane by squares and octagons. Here, 
all lines represent a resistor with resistance R. 

explicit expression for the resistance between two arbitrary nodes of any lattice structure 
which is a periodic tiling of d dimensional space with electrical resistors. We present 
several examples of nontrivial resistor networks that, to our knowledge, have not been 
studied in the hterature. Our Green's function method enables us to derive, for the 
first time, resistance formulas for the Kagome and dice lattices (see sections 13.31 and 
13. 4p . and for the decorated lattice (see section 13. 5p in terms of the resistances on the 
well studied triangular and square lattices, respectively. Our general method is not 
restricted to lattice structures made from identical resistors as is demonstrated in a 
simple, one dimensional example. Moreover, our method can be extended to a uniform 
tiling of a surface of a toroid or cylinder. The theory of perturbed lattices developed 
in reference [10] can easily be extended to such resistor networks as we discuss in this 
work. 

In condensed matter physics, often the tight binding model (see, eg, [71 Eni ED ESI 
|23lElll25]) is a very good approximation for calculating the electronic band structure or 
transport properties of crystalline matter, like the recently discovered graphene [26t [27]. 
In such an investigation, knowledge of the Green's function corresponding to the 
Schrodinger equation is essential. Our results may be used to calculate the wave 
functions at the lattice points for complicated lattice structures. 

The text is organized as follows. In section |2] we outlined the general formalism to 
calculate the two point resistance on an arbitrary periodic lattice. Several, non trivial 
examples are presented in section |3] together with a few analytical results regarding the 
resistance between nearby lattice points. The conclusions are drawn in section HI In the 
Appendices we listed the resistance formulas that provide mapping of the Kagome and 
dice lattice to the triangular lattice, and the decorated lattice to the square lattice. 
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2. General formalism 

Here we derive a general integral expression for the resistance between two arbitrary 
lattice points in an infinite regular lattice of any resitor network. 

First, we specify the d-dimensional regular lattice by its unit cell as it is common 
in solid state physics. The lattice point is given by r = niai + n2a2 + ■ ■ ■ + naaa, where 
ai, a2, . . . , a^; are the unit cell vectors in the d-dimensional space and ni, n2, . . . , n^^ are 
arbitrary integers. We assume that in each unit cell there are p lattice points labeled 
by a = 1,2, ... ,p. Now we denote any lattice point by {r, a}, where r and a specify 
the unit cell and the lattice point in the given unit cell, respectively. Figure [H shows 
an example for the lattice structure of a resistor network. We assume that all lines 
represent a resistor with resistance R (although in the general formalism this restriction 
has not been used). Figure [2] shows one possible choice for the unit cell of the lattice 
shown in figure [TJ There are four types of lattice point in each unit cell denoted by 
a = 1,2,3,4 in figure [2 




a, 



Figure 2. The unit cell (dashed lines) with unit vector ai and a2 of the lattice shown 
in figure[T] There are four lattice points in each unit cell labeled by a = 1, 2, 3, 4. 

Denote the electric potential at the lattice point {r, a} by Vo(r) and the net current 
fiowing into the network at the lattice point {r,a} by /a(r) (as one can see below it 
is convenient to treat the problem by assuming that the current enter at {r, a} from 
a source outside the lattice). Owing to Ohm's and Kirchhoff's laws the current /^(r) 
relates to the potential V^(r) by the following equations: 

Y,La^i^-^')V^ir') = -Ur), (2.1) 
r',/3 

where Lap{r — r') is a p by p matrix called the connectivity matrix or Laplacian matrix. 
For concrete examples see section [31 The minus sign has been introduced only for 
convenience. In fact, if there is only one lattice point in each unit cell, ie, p = 1, then 
La^ir — r') is the lattice Laplacian corresponding to the finite-difference representation 
of the Laplace operator used in the literature (see, eg, references [281 IS])- From the 
translational symmetry of the lattice structure it follows that the Laplacian matrix 
depends only on the difference r — r'. Equation (12. ip can be written in a more compact 
form using matrix notation: ^j./ L(r — r') V(r') = — I(r). Owing to the nature of the 
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connectivity between lattice points, the Laplacian matrix satisfies the following relation: 
L(r) = L"'"(— r), where T denotes the transpose of a matrix. 

We now take periodic boundary conditions and consider a lattice with 
A^i, N2, . . . , Nd unit cells along each unit cell vector ai, a2, . . . , a^. Thus the total number 
of unit cells in the d dimensional lattice is iV = Y[i=i the total number of lattice 

points is Np. Now, it is useful to introduce the reciprocal lattice. Owing to the periodic 
boundary conditions, the wave vector k in the reciprocal lattice is limited to the first 
Brillouin zone (BZ) and is given by 

where mi, m2, ■ ■ ■ are integers such that —Ni/2 < rui < Ni/2 for i = 1,2, ■ ■ ■ ,d, and 
hj are the reciprocal lattice vectors defined by ajbj = 27r5jj, i,j = 1,2, - ■ ■ ,d. Here we 
assumed that each Ni is an even integer, which will be irrelevant in the limit iVj — )■ 00. 
The mathematical description of the crystal lattice and the concept of the Brillouin zone 
can be found in many books on solid state physics [Sni EH ESI ESI EH 125] . 
We take the discrete Fourier transform of the current I(r): 

I(k)= 5^I(r)e-^ (2.3a) 

r 

and an analogous expression is valid for the potential V^(k). Then the inverse Fourier 
transforms are given by 

IW = ^E We^^ (2.36) 

keBZ 

and the same is valid for the potential V^(r). This expression can easily be proved using 
the well-known relation J2r e''"^''"^'^ = N5k,k' [201 ED 1221 E31 El ES] ■ 

Using the Fourier transform of the current and the potential distributions 
equation (12. ip can be rewritten as 

L(k)V(k) = -I(k), (2.4a) 

where 

L(k) = 5^L(r)e-'^^ (2.46) 
r 

Since L(r) = L"'"(— r) it is obvious that L(k) is a Hermitian matrix, ie, L(k) = L'*'(k), 
where + denotes the conjugate transpose of a matrix. 

For a given I(k) the Fourier transform V(k) of the potential distribution can easily 
be found from (I2.4ap : 

V(k) = G(k)I(k), (2.5a) 

where the lattice Green's function G(k) (in k-space) is given by 

G(k) = -L-^(k). (2.56) 

Since L(k) is a Hermitian matrix, the Green's function G(k) is also a Hermitian matrix. 
If there is only one lattice point in each unit cell, ie, p = 1, the lattice Green's function 
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G(k) is a one by one matrix and it is the lattice Green's function of the finite-difference 
representation of the Laplace operator used in the literature (see, eg, reference [5]). 

To measure the resistance between lattice points {ri,ai} and {r2,a2}, we assume 
that a current / enters at lattice point {ri,ai} and exits at lattice point {r2,a2}, and 
at any other lattice points the currents are zero. Therefore, the current distribution in 
this case can be written as 

^uir) = I {6i,^ai (^r.n - Su,a2 (^r.ra) , (2.6a) 

and its Fourier transform becomes 

r 

Using equations fl2.66l) and (12.5 ap we find 

V,ik) = 5^G^,/,(k) = / [G^,,(k)e-ii'^^ - G^,,(k)e-*^^] , (2.7a) 

V 

and the potential distribution in r-space becomes 

^^(r) = ^ E ^'^(k)e^^r = ^Y. [G',.,(k)e^'^(--) -G,.,(k)e^^(— )] . (2.76) 

keBZ keBZ 

Now the resistance between lattice points {ri,ai} and {r2,a2} is given by 

fl„,„,(,,r,) = l^«M^iyi5). (2.8) 
Then using equation fl2.76l) we find 

-^0102(1^1)^2) = E I^G'aiai(k) +G'a2a2(k) 

keBZ 

CrQ,-^Q,2 (k) e (jQ,2ai(k)e . (2-9) 

If we take the limit A^j — )■ 00 for all i = 1, 2, . . . , d then the discrete summation over 
k can be substituted by an integral [201 [211 [22| [231 [2^ [25] : 

1 /■ d'^k , , 

(2.10) 

ksBZ "^keBZ (27r) 
where Vq is the volume of the unit cell. Thus the resistance -Ra^a2(ri, r2) becomes 

f d'^k r 

-Raia2(ri) I"2) = ^^0 / G'cjiai(k) + 6^0202 C^) 

JkeBZ (27r) L 

-G',,,2(k)e-*^'"'~'"^^-G',2ai(k)eii^('-^-^^)]. (2.11) 

This is the central result of this work. Note that the resistance depends only on 
the difference r2 — ri which is a consequence of the translational symmetry of the 
resistor lattice. Since G(k) is a Hermitian matrix, it is clear that the first two and 
the sum of the last two terms in the above expression are real, thus the resistance is 
real, as it should be. Moreover, from equation (12. lip the following symmetry relation 

-Roiaall"!) 1"2) = -Ra2,ai (r2, Tl) IS alsO obvioUS. 
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It is clear that deforming the lattice structure does not change the resistance 
between two arbitrary lattice points if the topology of the lattice structure is preserved. 
The same deformation was used by Atkinson and Steenwijk [3| for the triangular 
lattice. Thus, one can always deform the lattice structure of any resistor network 
into a d dimensional hypercube in which all the unit cell vectors have the same 
magnitude and they are perpendicular to each other. This topologically equivalent 
d dimensional hypercube lattice is more suitable for evaluating in equation (12. lip the 
necessary integrals over the Brillouin zone since the Brillouin zone also becomes a d 
dimensional hypercube. Indeed, if we write r2 — ri in terms of the unit cell vectors: 
r2 — Fi = niai + n2a2 + ■ ■ ■ + n^Sid then the exponentials in equation (12. lip and the 
Green's function will depend on the variables Xi = kai, X2 = ka2, . . . ,Xd = ka^^, ie, 

G{xi, . . . ,Xd) = G(kai xi, . . . ,kad Xd). (2.12) 

Moreover, the integration variables k = {ki, k2, ■ ■ ■ , kd) in (12.1 ip can also be transformed 
to the variables xi, X2, . . . , x^ for which the limits of the integration are ±7r. The 
Jacobian corresponding to the transformations of variables cancels the volume Vq of 
the unit cell in (12. lip . Thus, the general result (12. lip for the resistance can be rewritten 
as a. d dimensional integral that is a more suitable form for explicit calculations: 





r dxd 


27r . 


Ln 2vr 


—i{nix\+n2X2 + 





'~'a\a2\-'^li ■ ■ • 1 ■'^dj ^ a2ai\-'^li • • • 1 ■'^dj ■ y^'^.^uj 

For any infinite regular resistor network the above result is an integral representation 
of the resistance between arbitrary lattice points and provides a practical way for 
determining the resistance in question. 

3. Applications of the general formalism 

In this section we demonstrate how this result can be applied to different lattice 
structures. We present examples of one, two and three dimensional cases as well. 

3.1. One dimensional resistor structure 

Perhaps, the simplest periodic tiling is the one dimensional lattice structure in which 
each unit cell contains two non-equivalent resistors as shown in figure |3l This lattice 

R, 7?, 7?, i?T 
-• i— •■ — • — !— • — -• — 



Figure 3. A periodic tiling in one dimension with two resistors Ri and i?2- In this 
lattice each unit cell contains two non-equivalent lattice points labeled by a = 1, 2, ie, 
p — 2. The length of the unit cell is a. 



structure is also an example of the case in which the connectivity between lattice points 
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in each unit cell is not necessary the same. As we mentioned before, our general 
formalism outlined above can be applied to such cases as well. Indeed, now there 
are two different resistances in the unit cell, Ri and i?2- 

In one dimension we can characterize lattice points by {r,a}, where r = ma and 
a = 1,2 (here m is an integer). Ohm's and Kirchhoff's laws for the currents /a(r) at 
site {r, a} (with a = 1,2) can be written as 

Vi{r)-V2ir) Viir) - V2{r - a) 



hir) 



h(r) 



Ri 



R2 



V^{r)-Vi{r)_^V2{r)-Vi{r + a) 



Ri R2 
Hence, the Laplacian matrix can be written as 



L(r) 



<5r,0 

Ri 


R2 


Sr.O 

Ri 


Ih 


Sr,Q _|_ 

Ri ^ 


<5r,-a 
R2 


Ri 


R2 . 



(3.1a) 
{3.1b) 

(3.2a) 



and its Fourier transformation reads 
L{k) 



1 

Ri 


1 

i?2 


i + 


g-ifco - 

R2 


. i + 


gifca 
R2 


1 

Ri 


1 

R2 . 



{3.2b) 



G{k) 



4 sin^ Y 



(3.3) 



Inverting the matrix L(A;) the Green's function defined by equation ( 12.561) becomes 

-Ri + R2 R2 + Rie~^^"' 

R2 + RiQ'^^"" Ri + R2 

Now, changing the variable k to x, the Green's function in equation fl2.13p (for d = 1) 
transforms to G{x) = G{ka — )■ x). 

The resistance between lattice points {0,a} and {ma, /3} can be obtained from 
equation (12.131) for the one- dimensional case ( d = 1): 

Rn{m) = R22{m) = {Ri + R2) f{m), (3.4a) 



Ri2{m)=RJ{m+l) + R2f{m), 

R2i{m) = Rif{m — 1) + R2f{m), where 

_ , , r dxl — cos mx 

f{m) 



(3.46) 
(3.4c) 

{3.4d) 

2n 1 — cos X 

The function f{m) can be evaluated using the method of residues and we find f{m) = 
\m\, if m is an integer. Thus, -Rii(m) = (i?i + R2) \m\, Ri2{m) = Ri\m + 1| + R2\m\ 
and R2i{m) = Ri\m — 1\ + R2\m\. These results can easily be understood from the fact 
that the current flows only between lattice point {0, a} and {ma, /?}, and between these 
points the resistances are in series. The two semi-infinite segments of the chain do not 
affect the resistance. 
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3.2. Tiling of the plane by squares and octagons 

Consider the resistor network shown in figure [H Hereafter, we assume that all lines 
represent a resistor with resistance R. Using Ohm's and Kirchhofi^'s laws the currents 
/Q,(r) at site {r, a} (with a = 1,2, 3, 4) can be written as 



h(r) 



Vo(r) 



VAr) 





R 


+ — — 


R 


+ — 




R 










;r)- 




(r - as) 




R 




R 




R 








-ai) ^ 


Vs{r) - V^ir) 




[r) 






R 




R 






R 








- V2ir + aa) 


V4 







R 



R 



R 



(3.5a) 
{3.5b) 
(3.5c) 
(3.50 



From these equations one can easily read out the Laplacian matrix: 



Lfr) 



1 

R 



— 3(5r,r 


Sr,0 


5r,ai 


Sr,0 


5r,0 


-3Sr,0 


^r,0 








-3(5r,0 


5r,0 


Sr,0 


<^r,-a2 


5r,0 


-36r, 



(3.6) 



Hence, the Fourier transform of the Laplacian matrix f l3.6p is given by 

-3 1 e-'^^ 1 
1 -3 1 e-'^'' 



L(k) = ^ 



1 -3 

^ika2 



(3.7) 



1 

1 e"^2 1 _3 

Inverting the above matrix L(k), the Green's function can easily be calculated from 
equation (12.5&I) and changing the variables we find 

G{xi,X2)= G(kai — )■ a;i, ka2 — )■ X2). (3.8) 

Since the resistance depends only on the difference r2 — ri we can take the lattice vector 
ri (at which the current / enters the network) to the origin. Specify the lattice vector 
r2 = fq = mai + na2 (at which the current / exits the network). Then using the main 
result 02.131) for d = 2, the resistance between lattice point {0,a} and {rg,/?} is given 
by 



Raii{m,n) 



^ dxi r dx2 r 



2ti 



27r 



Gaa{Xl,X2) + Gf)ii{Xi,X2) 



G«/3(xi, X2)e-^("^^i+"^^) - G^„(xi, X2)e'(™^^+"^^)J . (3.9) 

It can be shown that in the above expression one integral can be evaluated analytically 
by the method of residues [28] much in the same way as in reference [5] for the square 
and triangular lattices. 
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Applying this result we calculated the resistance for a few cases. For example, for 
the resistance between lattice point 1 and 2 that belong to the same unit cell, we obtain 

, -9 + 4cosa;i + 4cosa;2 + cos (xi - 0:2) /o-,..n 
f{Xi,X2) = — ■ . 3.10a 

— 7 + 3 cos Xi + 6 cos X2 + cos Xi COS X2 

Performing the integrations we find 

, , fl y2arctan (2V2)\ 
Ri2{0, 0) = I - + 1 (3-106) 

and numerically, i?i2(0,0) ^ 0.6385i?. 

Similarly, the resistance between lattice point 1 and 3 (when they are in the same 
unit cell) is: 

, , 3V2 ( 2arctan\/2\ 
i?i3(0, 0) = 1 — ^ j R ^ 0.8312i?. (3.10c) 

Owing to the symmetry of the lattice, -Ri4(0, 0) = -R23(0, 0) = -R34(0, 0) = -Ri2(0, 0), and 

i?24(0,0) = i?i3(0,0). 

In adjacent unit cells the two squares are connected by one resistor, and then one 
can ask what is the resistance between the two ends of this resistor. In our notation, 
this can be found by calculating, eg, i?3i(l,0). Using equation (13. 9 p we obtain 

fl3.(l,0)= (^l-^ + ^^iH£^^jB, (3.10<i) 

and numerically, i?3i(l,0) ~ 0.7229-R. For more examples, it is useful to apply the 
method of residues. 

3.3. The Kagome lattice 

The next example we consider is the so-called Kagome lattice structure [25] shown in 
figure m The unit cell of the Kagome lattice structure is shown in figure [5l Each unit 



Figure 4. The Kagome lattice structure of the resistor network. 

cell contains three lattice points, ie, a = 1,2,3 and the Laplacian matrix is a three by 
three matrix. In a similar way as in section 13. 2[ one can write down the equations for 
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Figure 5. The unit cell of the Kagome lattice with three lattice points, p = 3. 



the currents for each lattice point in one unit celL Then it is easy to show that the 
Fourier transform of the Laplacian matrix for the Kagome lattice is 

-4 l + e-'^' l + e-^'^^ 

^^gikai _4 ;^ ^ gik(ai-a2) _ (^^^^ 

I _|_ e''^2 I _|_ g-ik(ai-a2) _^ 

As we mentioned before equation f l2.12p the lattice can be deformed to a topologically 
equivalent one in which the unit cell is a hypercube, in this case, it is a square lattice 
(see figure [6]). Now, calculating the Green's function from equation (12.5 6p and changing 



L(k) = ^ 




Figure 6. The topologically equivalent lattice structure of the Kagome lattice. The 
unit cell is a square. 



variables according to equation fl3.8p one can use the general expression (13 .Qp to find the 
resistance between arbitrary lattice points on the Kagome resistor network. 

It is interesting to calculate the resistance between adjacent lattice points. From 
symmetry it follows that -Ri2(0,0) = -Ri3(0, 0) = i?23(0,0), and from equation (13.91) 
we find that the integrand of the sum of these resistances is a constant and equals 
3R/2. Therefore, in equation (13. 9p the double integral results in -Ri2(0, 0) + -Ri3(0, 0) + 
-^23(0,0) = 3i?/2 and the resistance between the nearest neighbor lattice points is 
-^12(0,0) = i?i3(0,0) = i?23(0,0) = R/2. Based on the superposition principle and the 
symmetry of the lattice structure the same result was found in references [291 130] . 

However, our general formalism enables us to find the resistance between arbitrary 
lattice points. Moreover, we show that the resistance i?Q,^(m,n) on the Kagome lattice 
can be expressed in terms of the resistances R^{m,n) on a triangular lattice. To this 
end, we recall the result for the triangular lattice as discussed in references [U |5] . The 
resistance between lattice points {0, 1} and {fq, 1} (here fq = mai + nSL2 and for the 
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triangular lattice there is only one lattice point in each unit cell which is labeled here 
by 1) is given by 

/"^ dxi dx2 
J -^fmn{xi^X2), where 

fU.„.,)= l-cos(.n..+n..) 

6 — cosa;i — cosa;2 — cos(xi — X2) 
The resistor problem for triangular lattice has been well studied in references [U E] 
including a few analytical results. Note that the above expression is slightly different 
from that in reference |S] since here the angle between the two unit cell vectors ai and 
Si2 is 60°, while in reference [5] it is 120°, but this different choice of the unit cell is 
irrelevant for the resistances. Now, starting from equation (13. 9p one can show that the 
resistance Rai3{m, n) on the Kagome lattice can be expressed in terms of the resistances 



R^{m,n) on the triangular lattice. The explicit expressions are given in Appendix A 

Such an interesting map between the Kagome and the triangular resistor networks 
can be understood qualitatively by applying the so-called triangular-star transformation 
to the Kagome lattice just as it was used by Atkinson and Steenwijk for honeycomb 
lattices [1]. However, the Green's function method provides a more systematic way to 
deduce the resistance formulas as it was shown for the honeycomb lattice in reference [5] . 
Indeed, mathematically the map between the Kagome and the triangular resistor 
networks is based on the fact that the determinant of the Laplacian (13. lip for the 
Kagome lattice differs from the denominator of the function /^„(xi,X2) only by a 
constant factor (as it was the case for honeycomb lattice [5]). This mathematical fact 
should have a topological origin, which can be a research topic in the future. 



Using these results (see Appendix A ) it is easy to find the resistance between the 
opposite corners of a hexagon on the Kagome lattice, eg, by calculating i?33(l,0) or 
i?22(0, 1) and we obtain R^^il, 0) = i?22(0, 1) = + ^) R ^ 0.8120i?. 

3.4- The dice lattice 

The dice lattice (figure [7]) is a periodic tiling of the plane by rhombi having 60° and 
120° interior angles and all vertices have degree 3 or 6. The dice lattice is a frequently 
used lattice structure in the literature. Recently antiferromagnetic behavior has been 
studied on the dice lattice [31]. One possible choice of the unit cell is shown in figure El 
it contains three lattice points. The dice lattice can be deformed to a topologically 
equivalent, two dimensional square lattice as shown in figure [9l 

Using the method outlined previously one can easily obtain the Fourier transform 
of the Laplacian matrix for dice lattice and it is given by 

-6 A* B* 



A -3 
B -3 



(3.13) 
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Figure 7. The dice lattice structure. 




Figure 8. Tlic unit cell of the dice lattice with three lattice points, p — 3. 




Figure 9. The lattice structure that is topologically equivalent to the dice lattice 
structure. 



where A = l + (^^' + e^'^^ ^nd B = (^^' + e''^^ ^ gik(ai+a2)_ 

Now, one can again use equation fl3.9p to calculate the resistance between arbitrary 
lattice points on a dice lattice. Much in the same way as in the case of Kagome lattice 
from equation (13 .Qp we find that the resistance between arbitrary lattice points on a dice 
lattice can again be expressed in terms of the corresponding resistance on the triangular 



lattice. The explicit expressions are given in [Appendix B[ The mathematical reason for 
this mapping is the same as in the case of Kagome lattice mentioned in section 13.31 

Here we present a few analytical results. For instance, the resistance between 
lattice point 1 and 2 (which are in the same unit cell) may be found by calculating 
-Ri2(0, 0) and we obtain a simple result i?i2(0,0) = R/2. Moreover, the resistance 
between the two ends of the shorter and longer diagonal of a rhombus on a dice 
lattice may be obtained by calculating -R23(0, 0) and -Rii(l, 0), respectively and we find 
i?23(0, 0) = (l + iJ) ^ 0.7393i? and Ru{l, 0) = R/2. 
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3.5. The decorated lattice 

Here we consider the so-called decorated which is a square lattice with extra lattice 
points in the middle of each bond as shown in figure [TH] together with a natural choice 
for the unit cell vectors and lattice points labeled by a = 1, 2, 3 in each unit cell. This 
lattice was used to study, eg, the phase diagram of a decorated Ising system [32] and 
recently the Hubbard model [33]. 









>3 

2 







1 a, 



Figure 10. The decorated lattice structure. Between each pair of full circles the 
bond represents a resistor with resistance R. The unit cell contains three lattice points 
labeled by a = 1, 2, 3, ie, p — i. 



In much the same way as in our earlier examples one can easily find the Fourier 
transform of the Laplacian matrix for the decorated lattice: 



Lfk) 



1 
R 



(3.14) 



-4 1 + e-'^^ 1 + e-'^^ 

1 + e'^^ -2 

1 + -2 

Then, the resistance between arbitrary lattice points can be obtained from 
equation (13.91) . However, again just as in the case of the Kagome and dice lattices here 
there is a close relation between the decorated lattice and the square lattice regarding 



the resistances. We derived explicit expressions (see Appendix CD for the resistance 
Rai3{jn,n) on a decorated lattice in terms of the resistances R^{m,n) (between the 
origin and the lattice point mai + nsL2 on a square lattice) given by 



R^{m,n) =R 



2n 



2lT 



■/mn(^ 1)^2), where 



cos(mxi + nx2) 



(3.15) 

2 — cosxi — C0SX2 

Regarding the resistor problem on square lattice see, eg, references [H [5] in which a 
few analytical results |1] and recurrence formulas [5] are presented. The mathematical 
reason for the mapping of the decorated lattice to the square lattice of resistor network 
is again based on the fact that the determinant det(L) of the Laplacian f l3.14p differs 
from the denominator of the function fmni^iy^^) only by a constant factor, similarly as 
in the case of Kagome lattice mentioned in section 13.31 
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Note that the first equation in Appendix C can easily be understood by a simple 
argument. To find the resistance Ru{m,n) one can disregard the middle points of 
each bond on the decorated lattice. Then the decorated lattice is equivalent to a square 
lattice in which the bonds have resistance 2R. Thus we find that i?ii(m, n) = 2R^{m, n) 
in agreement with that given by the first equation in Appendix C For example, 
i?n(l,l) = 4i?/7r. 



Using the results given in Appendix C we find that the resistance between lattice 
points 1 and 2 that belong to the same unit cell is i?i2(0,0) = 3i?/4. Similarly, 
the resistance between lattice points 2 and 3 (which are in the same unit cell) is 
i?23(0, 0) = (1 + l/vr) R ^ 1.3183i?. 



3. 6. Centered square lattice 



In sections 13.31 and 13.41 we showed that the Kagome and the dice lattice can be mapped 
to the triangular lattice, while in section [3751 it was shown that the decorated lattice is 
mapped to the square lattice. In this section we show another interesting mapping. In 
particularly, we find that the resistor network discussed in section 13.21 (see figure [1]) can 
be mapped to the so-called centered square lattice shown in figure [TTl 




Figure 11. The centered square lattice structure. Between each full circles the bond 
represents a resistor with resistance R. The unit cell contains two lattice points labeled 
by a = 1,2, ie, p = 2. 



Using Ohm's and Kirchhoff's laws we find easily that the Fourier transform of the 
Laplacian is 



L(k) 



R 



A B* 
B -4 



where 



A = — 8 + 2 (cos kai + cos ka2) , 
B = (l + e'*^^) (l + e^'^2). 
Then, the Green's function from equation ( I2.5fep becomes 



G(k) 



det(L) 



4 B* 
B -A 



(3.16a) 

(3.166) 
(3.16c) 

(3.17) 
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where the determinant of L(k) is given by 

i?det(L) = 28 - 12(coskai + coskas) 

— 4coskai coska2. (3.18) 

Note that this determinant is exactly the same as that of the Laplacian in equation ( 13. 7p . 
Therefore, the Green's function should have the same denominator in both cases. The 
numerator of the Green's functions in the two cases is different. However, the structure 
of the numerator of the Green's function for L(k) given in (13. 7p is the same as that 
for the centered square lattice. Thus, one can work out the resistance formulas (not 
presented here) analogous to those listed in the Appendices to find the exphcit mapping 
between the lattice structure shown in figures [1] and [TTl 

Using equation (13. 9 p we calculated the resistance analytically for a few cases: 

v/2arctan (72/2) 



i?ii(l,0) 



i?i2(0,0) 



TT 



V2 



arctan 



^22(1,0) 



-1 



1 9^2 



arctan (2^2 



R. 



(3.19a) 



(3.196) 



(3.19c) 



7T Ait 

and numerically, i?ii(l,0) ^ 0.2771/?, i?i2(0,0) ^ 0.3615/? and /?22(1,0) ^ 0.5651/?. 



3. 7. Tiling of the plane by squares and triangles 

To demonstrate how efficient our general formalism is for calculating the resistance in 
a resistor network, we present results for a more complex lattice. Consider a network 
which is a periodic tiling of the plane by squares and triangles shown in figure [T5J The 
unit cell may be chosen as can be seen in figure | 




Figure 12. The lattice is a periodic tiling of the plane by squares and triangles. 
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Figure 13. The unit cell of the lattice structure shown in figure [T^l There are eight 
lattice points in each unit cell, p = 8. 



The Fourier transform of the Laplacian matrix for the lattice shown in figure [12] 
can easily be found: 



L(k) 



1 

R 



-5 


A* 


1 


B* 


1 





B* 





A 


-5 


A 


C* 


1 





B* 





1 


A* 


-5 


1 





A* 





1 


B 


C 


1 


-5 





A* 





1 


1 


1 








-5 


1 


B* 


1 








A 


A 


1 


-5 


1 


1 


B 


B 








B 


1 


-5 


1 








1 


1 


1 


1 


1 


-5 



(3.20) 



where A 



ikai 



B 



e'^^ and C = e^^'^^-^i). 
It is difficult to find analytical result since the integrand in equation fl3.9p is a very 
complicated function of its variables. However, numerically we calculated the resistance 
-Ra/3(0, 0) (in units of R) between lattice points a = 1, . . . , 8 and (3 = 1, . . . ,S that belong 
to the same unit cell and find 



i?a/3(0,0) 





r4 



r2 rs r2 rg Vj 



r-r 

r2 ry 

rs rs ri 

r2 rs rg rg 

re rs rs rs r2 

rj r-j re rs r^ r2 



rs r2 rs re 

ri rs rs Tq r2 

re ^5 rs r2 

r2 r4 r2 

r2 ri 

r2 



(3.21) 



r3 r2 r2 r2 ri r2 

where ri ^ 0.3849, rs 0.4038, rg ^ 0.5108, r4 ^ 0.5396, rg ^ 0.5585, re ^ 0.5647, 
r7 ^ 0.6230, rg ^ 0.6631. It is interesting to note that the resistance Ru = and 
-Ri6 = ^6 are almost the same numerically. However, rigorously they are not the same 
since they are not related to each other by symmetry. 
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3.8. Body centered cubic lattice 

Finally, we present a non-trivial example for a three dimensional lattice of a resistor 
network. The simple cubic lattice has already been studied in reference [5] , and Glasser 
and Boersma calculated the exact values of the resistances for a few cases [31]. Consider 
a more complicated resistor network, namely the body centered cubic (bcc) lattice shown 
in figure [TH 









a. 











Figure 14. The body centered cubic lattice as a periodic tiling of three dimensional 
space. Besides the sides of the cube there are resistors between the center of the cube 
and its corners. There are two non-equivalent lattice points in the bbc lattice, one is 
at one of the corner of the cube and the other is at the center of the cube, ie, p = 2. 
The unit cell vectors are ai, 3.2 and a^. All lines represent a resistor with resistance 
R. 



As we demonstrated above, the Laplacian matrix can be obtained from Ohm's and 
Kirchhoff 's laws and we find that its Fourier transformation can be written as 



Lfk) 



1 
R 



A-U B* 
B -8 



, where 



A = 2 (cos kai + cos ka2 + cos kaa), 
B = (l + e"«^i)(l+e''^^)(l+e^'^^). 
Then, the Green's function from equation ( 12.5fep becomes 



G(k) 



1 



8 B* 
B U-A 



(3.22a) 

(3.226) 
(3.22c) 

(3.23) 



det(L) 

where the determinant of L(k) is given by 

i?det(L) = 112 — 16(coskai + coska2 + coskaa) 

- 8(1 + coskai)(l + coska2)(l + coskag). (3.24) 

The resistance between two arbitrary lattice points can be determined from 
equation fl2.13p for d = 3. For a few cases we find numerically that the resistances 
(in units of R) are /2i2(0,0,0) ^ 0.1945, i?ii(l,0,0) ^ 0.1481, i?ii(l,l,0) ^ 0.1651, 
i?ii(l,l,l) ^ 0.1717 and i?22(l,0,0) ^ 0.2657. It is interesting to note that the 
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resistance i?ii(l,0, 0) (resistance between the two ends of a side of the cube) is much 
less than that in the simple cubic lattice (in this case it is -R/3, see, eg, [5]). The 
physical reason for this difference is that comparing with the simple cubic lattice, the 
body centered cubic lattice provides more channels for the current flowing between the 
two lattice points. 

4. Conclusions 

In this work using the Green's function method we derived a general resistance formula 
for any infinite lattice structure of resistor networks that is a periodic tiling of space in 
all dimensions. Our general resistance formula was applied to several non-trivial resistor 
networks to demonstrate how versatile our approach is. For the Kagome and dice lattice 
we derived exphcit expressions for the resistances between two arbitrary lattice points 
in terms of the resistances on a triangular lattice of resistors. Similarly, we showed that 
there is a direct map between the decorated lattice of resistors and the square lattice 
of resistors. We pointed out that a mapping between the lattice structure shown in 
figures [1] and [11] can exist. We believe that such a mapping between different lattice 
structures of resistor networks has a topological explanation. In fact we think that there 
exists a classification of different resistor networks in terms of some classes. However, 
for a deeper understanding of this issue more work needs to be done. This problem 
could be a future challenge for physicists and mathematicians. 

Tiling of plane is common in the arts and its mathematical description based on 
group theory is well known in the literature. Under the title 'tiling' one can find 
numerous decorative and practical examples of possible tilings on the world-wide web. 
Replacing the lines by resistors in such tilings provides a wealth of examples for possible 
resistor networks not studied in the literature. Here, we presented examples that are 
well known and relatively simple to find analytical results for. However, our Green's 
function method is general and makes it possible to study very complicated lattice 
structures such as that discussed in section [3171 

Our work can be extended to study the classical lattice dynamics and the vibrational 
modes of atoms within the framework of the harmonic approximation. Similarly, 
description of the electron dynamics governed by the Schrodinger equation in the tight 
binding approximation could be another application of the Green's function approach 
outlined here in the case of resistor networks. 
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Appendix A. Relation between the Kagome and the triangular lattice of 
resistor networks 

In this appendix we list the results for the resistance Rai3{m,n) on Kagome lattice in 
terms of the resistances R^{m, n) on a triangular lattice: 



Rii{m, n 
Ri2{m,n 
Ri3{m,n 
R22{m,n 
R23{m,n 
i?33(m,n 



R 7 1 
- + - i?^(m, n) - - R^{m-l,n + l) + R^{m + l,n- 1) 
9 o D 



R 5 
"9^6 

R 5 
"9^6 

R 7 
~9^3 



R''{m,n) + R''{m + l,n) 



+ 6 



R''{m,n) + R^{m,n + 1) 
1 



1 

+ 6 



R 5 r 

"9^6 
R 7 



i?"" (m + 1 , n - 1) + i?"" (m, n + 1) j , 
R''{m + l,n) + R^{m-l,n + l)^, 
R^{m,n)-^ R''{m,n-1) + R^{m,n + 1)^, 

R^{m -l,n) + R^{m, n + 1)] , 



R^{m,n) + R^{m - + 1) 
1 r 



— + -R''{m,n) 



R'^im - l,n) + R'^im + l.n) 



9 3 ' ' ' 6 

The remaining resistances can be obtained from the symmetry relation Raj3{m, n) = 
Rf}a{-m, -n). 

Appendix B. Relation between the dice and the trianguleir lattice of 
resistor networks 

Here we list the results for the resistance Rap{m,n) on a dice lattice in terms of the 
resistances R^{m,n) on a triangular lattice: 



Rii{m, n 
Ri2{m,n 

Ri3{m,n 
R22{m,n 
R23{m,n 



-R^{m,n), 



R 1 
6+2 

R 1 
6+2 



R^{m,n) + R'^im + l,n) + R^{m,n + 1)^, 

R^{m + l,n) + R^{m,n + 1) + R^{m + l,n + 1) 



R 1 
3+6 



2i?^ (m, n) + 2i?^ (m + 1 , n) + 2i?^ (m, n + 1) 
+R''{m - l,n + 1) + R'^im + l,n - 1) + R'^im + l,n + 1) 
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R3s{m,n) = R22{m,n), 

where = 0, if m and n equal to zero, otherwise it equals to 1. The remaining 
resistances can be obtained from the symmetry relation Rap{m,n) — Rpa{—'fn-i —n). 



Appendix C. Relation between the decorated lattice and the square lattice 
of resistor networks 



Here we list the results for the resistance Ra/3{m,n) on a decorated lattice in terms of 
the resistances FP{m,n) on a square lattice: 

Rii{m,n) = 2i?°(m, n), 
R 

Ruim, n) = J + i?°(m, n) + R°{m + 1, n), 
R 

Risim, n) = — + R°{m, n) + i?°(m, n + 1), 

i?22(m, n) = - + 3 i?°(m, n) - - [i?°(m, n - 1) + i?°(m, n + 1) 
2 2 - 



R23{m,n) 



R 1 
2+2 



i?°(m, n) + i?°(m - 1, n) + R^{m - 1, n + 1) + i?°(m, n + 1) 

i? 1 r 

i?33(m,n) = - + 3i?°(m,n) - - i?°(m - 1, n) + i?°(m + 1, n) . 

The remaining resistances can be obtained from the symmetry relation Rap{m,n) 
Rl3a{-m, -n). 
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